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We revisit the analysis of the improved ladder Schwinger-Dyson (SD) equation for the dynamical 
chiral symmetry breaking in QCD with emphasizing the importance of the scale ambiguity. Previous 
calculation done so far naively used one-loop MS coupling in the improved ladder SD equation 
without examining the scale ambiguity. As a result, the calculated pion decay constant / T was less 
than a half of its experimental value f n = 92.4MeV once the QCD scale is fixed from the high 
energy coupling af s (A4z)- In order to settle the ambiguity in a proper manner, we adopt here 
in the present paper the next-to-leading-order effective coupling instead of a naive use of the MS 
coupling. The pion decay constant f n is then calculated from high energy QCD coupling strength 
ctf s (Mz) = 0.1172 ± 0.0020. Within the Higashijima-Miransky approximation, we obtain f„ = 85- 
106MeV depending on the value of ctf s (Mz) which agrees well with the experimentally observed 
value f w = 92.4MeV. The validity of the improved ladder SD equation is therefore ascertained more 
firmly than considered before. 
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I. INTRODUCTION 

The improved ladder approximation of the 
Schwinger-Dyson (SD) equation [jl], |j has been 
used widely for the analyses of the dynamical chi- 
ral symmetry breaking in QCD pi 0, 0, || and 
of models of dynamical electroweak symmetry 
breaking J?], ||. In the improved ladder approx- 
imation, the tree-level one-gluon-exchange dia- 
gram, with its QCD coupling strength being re- 
placed by the running MS one, is assumed to give 
a dominant contribution in the attractive force 
between quark and anti-quark (qq). 

The method actually succeeded in explain- 
ing many quantitative relations among the low 
energy hadronic data in QCD ||, ||, ||, [| U 
and established its qualitative validity. Recently, 
it has also been adopted to predict properties 
of QCD under extreme conditions (hot and/or 
dense QCD) @ f| H [jj, 0. 

It is known, however, that a naive use of 
this approximation, combined with the high en- 
ergy QCD coupling such as af s (Mz), predicts 
the pion decay constant significantly smaller 
than its experimentally observed value = 
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92.4MeV. §, [| (See, for example, our Figure [§ 
where we obtained f v of order 30MeV for the one- 
loop running of the MS coupling.) Note that the 
decay constant determines the scale of chiral 
phase transition in QCD. Does this discrepancy 
indicate the existence of substantial non-ladder 
contribution in the driving force of the dynamical 
chiral symmetry breaking? If this is so, the suc- 
cess of the improved ladder approximation should 
be regarded as an accidental coincidence, without 
understanding deeply the nature of the driving 
force of the dynamical chiral symmetry breaking. 

Failing to explain the experimental value of f n 
from the high energy QCD coupling strength, the 
improved ladder SD equation is often regarded as 
a viable phenomenological model of low energy 
QCD, in which we can freely tune its coupling 
strength so as to make consistent with the 
experimental value. We emphasize, however, in 
order to perform trustful analysis of QCD under 
extreme conditions, the low energy models need 
to be connected smoothly with the high energy 
QCD. The discrepancy of thus causes a seri- 
ous trouble in the analysis of the critical behavior 
of hot and/or dense QCD. 

In this paper, we point out that the previ- 
ous calculations so far made for overlooked 
the existence of scale ambiguity |15[j in the im- 
proved ladder SD equation. In order to settle 
the ambiguity, we introduce a concept of effec- 
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tive coupling [[l6[ [l7| in the analysis of the im- 
proved ladder SD equation of QCD. (The idea 
similar to our effective coupling was proposed in 
the analysis of SD equation by Ref.|l7j for dif- 
ferent purpose to ours. It was also adopted by 
Ref. jll| [E| for the analysis of SD equation in 
gauge theories with extra dimensions.) The effec- 
tive coupling is calculated using the background 
field method |^0| in this paper. We next calcu- 
late numerically the improved ladder SD equa- 
tion and obtain the value of the pion decay con- 
stant /jr using the high energy QCD coupling 
af s (Mz) as an input parameter of the analysis. 
Within Higashijima-Miransky approximation jjj, 
we obtain f„ = 85-106MeV depending on the 
value of af (M2) = 0.1172 ± 0.0020 |l). The 
agreement of the calculated with the experi- 
mentally observed value f v = 92.4MeV (or the 
value in the chiral limit m u — = m s = 0, 
/tt = 86MeV j22j) is rather impressing. The dis- 
crepancy of is resolved in our analysis. We 
thus ascertained the validity of the improved lad- 
der SD equation more firmly than considered be- 
fore. 

We organize the present paper as follows: In 
Section || we give a brief derivation of the SD 
equation and point out the problem of the scale 
ambiguity. In Section |n| the concept of effective 
coupling is introduced. The behavior of QCD 
coupling and its regularization are discussed in 
Section IV. In Section [v| we give our results 
of numerical analysis within the Higashijima- 
Miransky approximation to the angular integral 
of the SD equation. A numerical analysis without 
the Higashijima-Miransky approximation is per- 
formed in Section ^ using the non-local gauge 
fixing parameter method. Section VII is devoted 
for summary and discussions. 



II. LADDER SCHWINGER-DYSON 
EQUATION 

We outline here a derivation of the ladder SD 
equation based on Ref. J23[ |24| , |25| ]. We first in- 
tegrate out the gluon field from the QCD La- 
grangian at the tree level. We then obtain a bi- 
local interaction model, 



5bl = Jd 4 x ipifiip + — J d A x\d A X2~* 

x(^T>) Xl (^ 7 ^ a V0sA^l-*2),(l) 



with gluon propagator D^ v (x) being given by 
D^ix) = 



d 4 k 



•(2) 



Here g s and £ are the QCD coupling constant 
and the gauge parameter, respectively. In Eq.(Q) 
we have neglected six or higher quark interac- 
tion terms. The ladder SD equation for the dy- 
namical chiral symmetry breaking can be derived 
from Eq. (Q) (bi- local four fermion interaction) in 
a similar manner to the gap equation in the NJL 
model (local four fermion interaction). We ob- 
tain 



A{ V) = 1 + % 
A(-q 2 ) 



d 4 q 
(2tt)H 



-A 2 (-q 2 )q 2 + B 2 (-q 2 



(1 + 



p-q 
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B(-q 2 ) 
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(3b) 



with C F = T a T a = (N 2 - l)/(2N c ) = 4/3 be- 
ing the Casimir of the fundamental representa- 
tion of SU(N C = 3). The functions A(-p 2 ) and 
B(—p 2 ) in Eq.(||) are quark wave-function and 
mass-function, respectively. The quark propaga- 
tor 5* is given in terms of these functions, 



iS- 1 (p)=i>A(-p 2 )-B{-p 2 ). 



(4) 



Nonvanishing mass-function B ^ implies the 
dynamical chiral symmetry breaking in QCD. It 
is therefore regarded as an order parameter of 
this system. 

Performing the Wick rotation and the angular 
integrals, we obtain 



A(x) 



B(x) 



1 + 



A 2 

dy 



yA(y) 



C F / dy 



yA 2 (y) + B 2 (y) 
yB(y) 



K A (x,y), 



o y yA 2 (y) + B 2 (y) 



(5a) 

K B (x,y). (5b) 



Here we have introduced the ultraviolet (UV) 
cutoff A in the SD equation. It is known that 
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the result of the SD equation is insensitive to the 
cutoff A if we take sufficiently large A in QCD. 
The integral kernels Ka and Kb of the SD equa- 
tion Eq. (0) are given by 



K A (x,y) = / d9 sin 2 9 x 



xa. 



(3-0 



xy cos ( 



z 



K B {x,y) 



1 

2^ 



-2(1 -Of sin 2 



(6a) 
.(6b) 



with z and a s being the square of the gluon mo- 
mentum, 



y — 2^/xycos9 



and the QCD coupling 



El 
4tt' 



(7) 
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respectively. 

The ladder SD equation Eq. (0) combined with 
the kernel Eq.(^j) does not take account of the 
running of QCD coupling strength, however. In 
order to make the ladder SD equation consis- 
tent with the renormalization group (RG) consid- 
eration, the bi-local four-fermion interaction in 
Eq. (|TJ) needs to be modified to include its running 
effects. A widely adopted prescription for such a 
purpose is to replace the QCD coupling g s with 
the MS renormalized one g s (/z) where /i is identi- 
fied with the gluon momentum sfz. jl], ^6|, [2?], 

The prescription is not unique, however. Ac- 
tually, we can use g 2 {\i = c^fz) (c ^ 1, c> 0) 
instead of the conventional choice g 2 {n = \fz) 
(scale ambiguity Jl5|). Even if we adopt c ^ 1, 
the solution of the SD equation is shown to 
be consistent with the RG. Moreover, as we 
stated previously, naive analysis based on the im- 
proved ladder SD equation with c — 1 predicts 
fn much lower than its experimentally observed 
value / w = 92.4MeV. g § 

In order to calculate the pion decay constant 
from the high energy QCD coupling a s (M z ), 
we thus need to resolve the scale ambiguity in 
the context of the improved ladder SD equation. 
We consider this problem in the next section by 
introducing the concept of effective coupling. 



III. EFFECTIVE COUPLING IN QCD 

Before investigating the scale ambiguity of the 
improved ladder SD equation of QCD, we first 
consider simpler model, the strongly coupled 
QED with N flavor of massless fermions |^9|, [30| . 

In this model, ladder SD equation can be de- 
rived from a bi-local four-fermion interaction 

WWk (^) X2 D^ B (x 1 - x 2 ), (9) 

which is induced by the photon propagator 



d k a 

. e -»*-*^£ ( 9ll „ + (gauge fixing)) 



(2tt)< 



k 2 



(10) 



with g e being the QED coupling constant. 

Thanks to the Ward-Takahashi (WT) identi- 
ties of QED, it is enough to calculate vacuum 
polarization function (loop correction to the pho- 
ton propagator) for the evaluation of the bi-local 
four-fermion interaction Eq.fl9J). At the one- loop 



level, the g 2 /k part of Eq.(10) is then replaced 

as 

A 9e 
k 2 



(1 



g 2 nQ ED( fc 2)) fc 2 
1 



1 



(11) 



with g e yfs being the MS renormalized QED cou- 
pling. The MS renormalization is performed in 
the second line of Eq.(ll). It is straightforward 
to calculate the MS regularized vacuum polariza- 
tion function Il2S- D : 



MS 



n S D (^ 



/0 = ~ 3 N 



In- 



(12) 



™° 6 (J,* 6 

We compare Eq. ( p"l| ) with the prescription of the 
naive improved ladder approximation, 



El 

k 2 



(13) 



The naively improved formula Eq.(fl3|) differs 
from our one-loop formula Eq. (|ll"| ) due to the ex- 
istence of nonvanishing finite part in the vacuum 
polarization, 

/tt^ 20 



U QED (k' 2 u 
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-N, 



-k A > 0. 



(14) 
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One method to resolve the scale ambiguity is to 
use the renormalization scale /i so as to minimize 
the finite corrections. In the present model, we 
find that the choice 



(15) 



eliminates the finite part 



n: 



'MS 



(fc 2 , M = v/|fc 2 |exp -- )=0. (16) 
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The other method is to use the effective cou- 
pling in the improved ladder SD equation. The 
effective coupling is defined by 



1 



1 



nQ| D (fc 2 , M ). (17) 



The ^-dependence in the RHS of Eq.(|l7|) can- 
cels at the leading-order. This method has an 
advantage compared with the former one that 
it can easily deal with massive particles in loop. 
We therefore adopt the effective coupling method 
hereafter. 

We now turn to the problem of the scale am- 
biguity of the improved ladder SD equation of 
QCD. As we have done in the case of QED, we 
need to calculate loop corrections to the bi-local 
four-fermion interaction Eq. (jl|) in QCD. A diffi- 
culty arises, however, in the case of QCD, where 
Slavnov- Taylor (ST) identities hold instead of the 
QED-like WT identities. Due to the lack of the 
QED-like WT identities, evaluation of the vac- 
uum polarization function is not enough to renor- 
malize the QCD coupling constant. 

Pinch technique (PT) method developed by 
Cornwall and Papavassiliou [EH], Q may be a 
hopeful candidate to resolve the scale ambigu- 
ity in the improved ladder SD equation of QCD. 
The effective coupling in the (S-matrix) PT is 
extracted from the on-shell qq scattering ampli- 
tude. Since the renormalization scale dependence 
should cancel in the on-shell amplitude, the scale 
ambiguity in the qq scattering amplitude is auto- 
matically resolved in the PT. In addition, it does 
not depend on the choice of the gauge parameter. 
We note, however, that these great features of the 
PT effective coupling is assured only when con- 
tributions of the pinch part of the ladder type 
diagram Figure ^ is included in the amplitude. 
The pinch part of Figure [l] gives non-zero am- 
plitude in gauges other than the Feynman gauge 



£ = 1. The ladder type diagrams are resummcd 
to all orders in the ladder SD equation. So, some 
portion of the ladder contribution will be doubly 
counted if we simply adopt the PT effective cou- 
pling in the analysis of the improved ladder SD 
equation with £ ^ 1. 



FIG. 1: Ladder-like diagram which is included in the 
definition of the PT effective coupling. 

We therefore adopt here a slightly different 
choice, the background gauge fixing method [pfj|| , 
where the QED-like WT identities hold even in 
the case of QCD. Thanks to the naive QED-like 
WT identities, the QCD coupling can be easily 
renormalized by the calculation of the vacuum 
polarization function in this method 
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1 



n^ D (fc 2 ,^), (is) 



l ' </D (fc 2 , fi) being the MS regularized vac 



with n 



MS 



uum polarization function. The MS renormalized 
QCD coupling is denoted by g s y[g- Using the one- 
loop finite part of n^2. D (/c 2 , fj,), it is easy to show 
that the ^-dependence cancels in the Eq.(|l8|) at 
the leading-order. We note that the coupling 
Eq.disl) can be free from the double counting 
problem in the ladder SD equation, since it does 
not include contributions from the ladder type 
diagrams. 

After a straightforward calculation, we obtain 



C G 
-2 



4J ( f(fe 2 , /U )-/f(A; 2 , M ) 

(l-£ bg ) 1 (1 -£ bg ) 2 
(4tt) 2 4 (4tt) 2 

Nf 

8T R Y,I?(k 2 ,m},fi). (19) 
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with C G = N c = 3, Tr = 1/2 for QCD 
(SU(N C — 3) gauge theory) and £b g being the 
gauge fixing parameter in the background gauge. 
The number of flavors is denoted by Np. Here 
the functions Iq 1 , 1^ are given by 



(4 7 r) 2 / ( f(fc 2 , / i) 
(47r) 2 /f(fc 2 ,/i) 



In- 



-k 2 
M 2 



2, (20a) 



1 — k 2 ^ 

— 4 < 20b > 



The function 1% represents the quark loop con- 
tribution. We obtain 



1 



(47r) 2 ^(fc 2 ,m 2 ) = --ln^ 



X 



tanlT 1 A, 



/ 



3 f 



— tanlT 1 A/ - 1 
x f 



with 



-k 2 



4m 2 



k 2 ' 



(21) 



(22) 



For a massless quark m/ = 0, the expression 
Eq.(|l]) reads 

(4.) 2 /«(fc 2 ,m 2 =0) = -iln-^ + 1 |. ( 2 3) 



Note here that the Eq.(18) depends on the choice 
of the gauge fixing parameter £b g - It is also 
known that the effective coupling Eq.(|l8|) with 
£bg = 1 coincides with the PT effective cou- 
pling. 
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Such a £bg-dependence cancels with the ladder- 
type contributions and other vertex corrections 
in the one-loop amplitude. The ladder-type dia- 
grams are resummed in the analysis of the ladder 
SD equation with a particular gauge parameter 
£. It is therefore plausible to take £bg = £■ As 
we will show in section [v|, the most convenient 
gauge parameter is £ = in the analysis of the 
ladder SD equation with Higashijima-Miransky 
approximation. We will therefore take £bg = 
in the analysis of following sections. We will also 
use £bg = 1 (the PT effective coupling) in order 
to make a comparison between them. 

We comment on the choice of the renormal- 
ization scale fi which minimizes the finite correc- 
tion Eq.([L9j). It is easy to find, for Np = 0, the 



finite correction vanishes with the renormaliza- 
tion scale n = a/|F| exp(-205/264) for £ bg = 
(fi = exp(-67/66) for £b g = 1). 



IV. BEHAVIOR OF RUNNING 
COUPLING 

We next evaluate numerically the behavior of 
the running couplings of QCD. Since we are in- 
terested in the dynamical chiral symmetry break- 
ing, we consider the running below the scale of 
b-quark mass nib = 4.3GeV. 

In our analysis, we use the PDG average J2T) 
of the high energy QCD coupling a^ s (fj, — Mz), 



af s (M z ) = 0.1172 ±0.0020, N F = 5, (24) 

where the MS coupling is defined in QCD with 
Np — 5 flavor of quarks. Eq.(|3) corresponds to 



af S (m b ) = 0.2197 ±0.0075, N F = 4, (25) 

at the next-to- next-to-leading-order (NNLO), 
and 

a^(m 6 ) = 0.2188 ± 0.0074, N F = 4, (26) 

at the next-to- leading-order (NLO) approxima- 
tion. We note that the difference between NLO 
and NNLO is of negligible order at n = m&. 

The renormalization group equation (RGE) 
of MS coupling is solved numerically below the 
scale of nib using Eq.(p5[) as its boundary con- 
dition. As we noted before, the difference be- 
tween Eq.(fj5|) and Eq.(p6|) is insignificant. We 
assume m u = nid = m s = for three light 
quarks. The Np = 4 MS RGE is adopted for 
/i > m c — 1.3GeV, while Np = 3 is used for 
[i < m c . Figure || shows the running behavior 
of one- and two-loop MS couplings. Note here 
that the two-loop effects become important for 
[i < lGeV. 

We emphasize here that we do not use Eq. ( |24| ) 
directly as a boundary condition of the one-loop 
RGE. Actually, the two- loop RGE effect between 
Mz and m& is of non-negligible order. If we adopt 
directly Eq. ( |24| ) as the RGE boundary condition, 
the one-loop RGE leads to the value of Q.y s ( TOb ) 
sizably smaller than the values of Eq.(P5[) and 
Eq.(|2q). The behavior of one- loop MS coupling 
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FIG. 2: Running of the MS (one- and two-loop) 
and effective (£b g = and £ bg = 1) couplings in 
QCD. The boundary condition of RGE is taken 



to be af b (mi,) 



0.2197, which corresponds to 



af s {M z ) = 0.1172 at the scale of M z . 



with this boundary condition is then quite differ- 
ent from the result of Figure ^. 

We next investigate the behavior of the effec- 
tive coupling. By using Eqs.(18- 2l]j2^ ) it is easy 
to show 



4tt 
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Here Iliight and Il c harm are defined as 
Iliight = — Gq 



11 -fc 2 
— In — 
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3 T ^ 3 



In 



-fc 2 



(28a) 



— tanh 1 X c - 1 



1 

X 2 



X c 3 c 



(28b) 



with JVg being the number of light flavors Ng — 3. 
The parameter X c is defined as 



Am 2 



(29) 



The MS coupling af s ( r i) in Eq.@ is the cou- 
pling at the scale /i > m c and thus denned in 4 
flavor QCD. For a^ifi < m c ) (3 flavor QCD), 
the In m 2 //i 2 term is just missing, 



4tt 



af (-fc 2 , fx < m c ) 



4ir 



vMS 



n 



light n c h arm 



(30) 



By using the expansion of tanh 1 X c around 



tanh X r 



X r 



X c 



1 



-xi 



:Xi 



,(31) 



it is easy to check the decoupling of the c-quark 
effect in the m c — > oo limit in Eq.([30|) 

We show the behavior of the effective coupling 
for £b g = and £b g = 1 in Figure [|. Since we 
want to calculate them at NLO level, we adopt 
the two- loop MS coupling in Eq. @ and Eq. (|§) . 
We also assumed fi = y/—k 2 in the plot. We 
note that the effective coupling a^f is signifi- 
cantly larger than the MS one. 



CO 



a 



a i 



0.5 0.6 0.7 0.8 0.9 



z [GeV] 



FIG. 3: Typical behavior of the IR regularized effec- 
tive coupling. 



Both the MS coupling af s (p = V^W = yfz) 
and the effective coupling af (-fc 2 = z, n = 
\J—k 2 = s/z) diverge in the infrared (IR) re- 
gion. We thus need to regularize them in order to 
solve the improved ladder SD equation with these 
coupling strengths. The IR regularized coupling 
a s (z) is regularized such that da s (z)/dz is con- 
tinuous for A 2 > z > 0. (A is the UV cutoff of 
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the SD equation.) In this paper, we adopt an IR 
regularization scheme || , 

a s (z = A 2 e t ) 

oto for t < to, 

a -^C{t~t ) 2 for*o<*<*i,(32) 
af(z = AV) for t x < t, 



with tt, to, C being determined by 
af (z = AV 1 ) = a x , 

to = h + 



2(a - 


- ai) 


d eff 




z— of 1 




az 


z=A 2 e*i 



c = 



1 



2(a - ai) 



z Tz a ° 



1 2 



z=A 2 e*l 



(33a) 
(33b) 

(33c) 



We also adopt the same IR regularization for the 
MS coupling. The typical behavior of the IR reg- 
ularized coupling a s is depicted in Figure |. Note 
here that our IR regularization Eq.(p2[) contains 
two IR coupling parameters ao and a.\. We re- 
gard the result of the SD equation reliable only 
when they are not too sensitive to these IR pa- 
rameters. 



V. CALCULATING WITH 
HIGASHIJIMA-MIRANSKY 
APPROXIMATION 

The angular integral Eq.(|^) cannot be per- 
formed in an analytical manner. In this section, 
we consider so called Higashijima-Miransky M 
approximation, in which a s (z) in Eq.(^) is re- 
placed as 



,(z) -» a s (max(aj,y)). 



(34) 



We can then analytically perform the angular in- 
tegral. We obtain 



£ r y 

K A (x,y) = — a s (x)-6(x-y) + (x<-+y) 

47T L X 



K B {x,y) 



3 + g 

47T 



a s (x) 



(35a) 

(x - y) + (x <-> y) . 

(35b) 



It is known that the existence of non- vanishing 
correction to the wave function A(x) makes 
the ladder approximation inconsistent with the 



gauge symmetry. We thus take Landau gauge 
£ = here in our numerical calculation of the 
improved ladder SD equation. 

The quark mass function B(x) can be regarded 
as an order parameter of the chiral phase transi- 
tion in QCD. Non-zero solution of the SD equa- 
tion B{x) thus indicates dynamical chiral sym- 
metry breaking and implies the appearance of the 
Nambu-Goldstone (NG) boson, the pion. The 
decay constant of pion / w can be related with 
the mass function B{x). 

In the following numerical analysis, we use the 
Pagels-Stokar formula for the pion decay con- 
stant [p5|, 



,2 _ 

JPS — A_r> 



Nc 

4tt 2 



dxx Adx a . (36) 



(x + B 2 (x)Y 



The Pagels-Stokar formula Eq.(J36[), however, 
does not fully take account of the pion wave- 
function of the ladder approximation. The for- 
mula is therefore not ladder-exact one. We thus 
relate it with the ladder-exact /„• by 



I-k — Afps/ps- 



(37) 



Fortunately, the Pagels-Stokar formula is known 
to work extremely well within the Higashijima- 
Miransky approximation of the ladder SD equa- 
tion. [p| We thus assume 



A^ps ~ 1, 



(38) 



in this section. 

We are now ready to perform numerical anal- 
yses of the improved ladder SD equation. The 
UV cutoff A is taken to be A = m&. We adopt an 
algorithm similar to the one described in Ref.| 18 
in our numerical analysis. 

Let us first confirm previous results of Ref. 0, 
Q, where the momentum dependence of a s (z) 
is given by the MS coupling with the renormal- 
ization scale /i being naively identified with y/z. 
Figure |]a shows the result of f w for the one- and 
two-loop MS coupling. We note that the cal- 
culated fn is extremely stable over wide range 
of the IR coupling regularization parameters ao 
and a\. The result is also consistent with pre- 
vious one that the calculated f v is significantly 
smaller than its experimentally observed value 
U = 92.4MeV. 

As we expect from the running behavior of a s 
(Figure 0) , the f„ calculated from two- loop MS 
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FIG. 4: Results of (a) the pion decay constant f n and (b) the "constituent quark mass" m const with use of 
the one- and two-loop MS couplings. The Higashijima-Miransky approximation is adopted in the kernel of the 
ladder SD equation. The boundary condition of RGE is assumed to be af s (mj) = 0.2197, which corresponds 
to af s (M z ) = 0.1172 at the scale of M z . iVps is taken to be 1. 



a s is larger than the value from the one-loop MS 
coupling. Although three- or higher-loop running 
effects may further enhance the value of /„, we 
expect that they are not enough to reproduce the 
value fa = 92.4MeV, since the higher-loop MS 
coupling is still well below the effective coupling 
at the scale ~ 700MeV, where the dynamical chi- 
ral symmetry breaking takes place in the analysis 
of the ladder SD equation. || 

We also emphasize that the mass function at 
the zero momentum B(0) depends significantly 
on the choice of the IR parameters ao and a.\. 
We find, however, that the "constituent quark 
mass" defined as || 

mconst = B(4m^ onst ), (39) 

is rather insensitive to aa and a\. (See Fig- 
ure ||b.) Our result on m const is also consistent 
with the previous analysis. 

We next consider the improved ladder SD 
equation combined with the effective coupling 
Eq.@ and Eq.@. We adopt two-loop MS 

coupling as aip(^) in Eq.@ and Eq.@. The 
leading-log dependence on fi cancels in the def- 
inition of cVg as we noted before. We take 
af(- k 2 — z,fi = i/z) in order to take account 
of the next-to-leading-log effects. 

We show our results on / w and m const in Fig- 
ure ||. Again, the results are quite insensitive to 
the variation of IR regularization parameters ctQ 
and a\. We also note that the calculated / w from 
a® with £bg = agrees very well with its exper- 
imentally observed value / ff = 92.4MeV. The 
result suggests the validity of our method to set- 



tle the scale ambiguity in the improved ladder 
SD equation. It may also imply importance of 
the ladder diagrams in the mechanism of the dy- 
namical chiral symmetry breaking of QCD. 

If we adopt £bg = 1 in the effective coupling, 
on the other hand, the calculated f„ becomes 
sizably larger than 92.4MeV (of order 120MeV). 
This deviation of /„■ in £b g = 1 may be under- 
stood as an indication of the double counting of 
the ladder diagrams as we described before. This 
point should be investigated further in future. 

In order to figure out the effects of the two- loop 
RGE in the effective coupling, we next perform 
analyses using one-loop cv^ s (/i) in Eq.(p7|) and 
Eq.(|30|). In this case, the //-dependence is com- 
pletely canceled in the effective coupling. Our 
results are shown in Figure || These results are 
sizably smaller than the results of Figure ||. The 
effect of two-loop RGE is therefore sizable in the 
effective coupling. Although the f % at £b g = 1 
agrees well with its experimental value 92.4MeV, 
the agreement should be regarded as an acciden- 
tal coincidence. Actually, the results shown in 
Figure ^ sizably decrease if we adopt Eq. ( p4^ ) di- 
rectly, instead of Eq.(p5|), as the boundary con- 
dition of the one- loop RGE. 

We have so far used the central value of 
Eq.(||), af s (m b ) = 0.2197 as an input param- 
eter of the QCD coupling. Larger the af s (m(,), 
larger f n is obtained in our framework, however. 
The uncertainty of the high energy QCD cou- 
pling thus leads to an uncertainty of our results 
on / w . There also exists uncertainty coming from 
the choice of the IR regularization parameters qq 
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FIG. 5: Results of (a) the pion decay constant and (b) the "constituent quark mass" m const with use 
of the effective couplings (two- loop RGE + finite part). The cases with £b g = and £b g = 1 are shown. 
The Higashijima-Miransky approximation is adopted in the kernel of the ladder SD equation. The boundary 
condition of RGE is assumed to be af s (mi,) = 0.2197, which corresponds to af s (Mz) = 0.1172 at the scale 
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FIG. 6: Results of (a) the pion decay constant f v and (b) the "constituent quark mass" m const with use of 
the effective couplings neglecting their two-loop RGE effects. The cases with £b g = and £b g = 1 are shown. 
The Higashijima-Miransky approximation is adopted in the kernel of the ladder SD equation. The boundary 
condition of RGE is assumed to be af s (mi,) = 0.2197, which corresponds to af s (Mz) = 0.1172 at the scale 
of Mz- Nps is taken to be 1. 



and a% . Since our results (Figures ||-||) are rather 
insensitive to «o & n d ol\, however, the bulk of un- 
certainty comes from the high energy QCD cou- 
pling here within the Higashijima-Miransky ap- 
proximation. 

We show the uncertainty of our results in Fig- 
ure 0. The high energy QCD coupling is taken in 
range of Eq. (p5|) . The IR regularization param- 
eters are also varied in range of 27r < ao < 4.7T, 
0.4 < a\/ao < 0.8 in the plot. The experimen- 
tally observed value of / ff = 92.4MeV and the 
decay constant in the chiral limit f v = 86MeV 
are also plotted in Figure |7| It is very impressing 
that the calculated value of / w agrees well with 
its experimental value for the effective coupling 



with £bg = 0. 



VI. CALCULATING f w WITH 
NON-LOCAL GAUGE 

We next try to perform the angular integral 
of the kernel Eq.(|fj) numerically without using 
the approximation Eq.(j34|). Unfortunately, the 
function A(x) receives non-vanishing correction 
even if we take £ = in this case. In order to keep 
A = 1, we need to use so called "non-local gauge" 
fixing in which the gauge parameter £ is regarded 
as a function of the gluon momentum. |56| [5^] 

It is known that A = 1 can be achieved when 
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we take a particular form p7| of function £, 

= 2 ^ / \ 

z 2 a s {z) 
We then obtain 



j dzz 2 -^-a s (z). (40) 



K A (x,y) = 0, 

K B (x,y) = [ desm 2 
27r Jo 



(41a) 
(3 + g(z))a a (s) 

z 

(41b) 



with z — x + y — 2-^/57/008 9. 

We first start with naive improved ladder case, 
in which a^ s (/i = %fz) is used as the running 
coupling of the SD equation. The naive im- 
proved ladder SD equation with one-loop a^ s 
has been investigated in full detail by Rcf.g. 
Unlike the case with Higashijima-Miransky ap- 
proximation, it was found that the Pagels-Stokar 
formula Eq.([56|) gives sizably larger f v than the 
ladder-exact value of in the non-local gauge 
fixing method. The factor ./Vps in Eq.(|37]) is then 
smaller than unity. We extract here the value 
of iVps from Ref.Q for various IR coupling pa- 
rameters ao, o.\. The results are summarized in 
table |. We find A^ps is of order 0.8 and it is a 
decreasing function of ao, ot\. We assume here 
^Vps = 0.8 irrespective to ao,«i throughout in 
this section. 

Figure |^ shows our results of the naive im- 
proved ladder SD equation with one- and two- 
loop MS couplings. Our result of the one-loop 



Uo/lT Ctl/n 


Aps 


2.3 0.74 


0.85 


3.1 0.89 


0.83 


4.4 1.1 


0.80 


7.2 1.5 


0.76 



TABLE I: The value of Aps for various ao and ai 
with one-loop MS RGE. 



MS coupling is consistent with the analysis of 
Ref.||. Although both f n and m const are rela- 
tively stable over the variation of ai/ao, their de- 
pendence on ao is of non-negligible order. We are 
thus not able to calculate / w in a reliable manner 
within the non-local gauge fixing method. The 
size of the calculated f v is, however, significantly 
smaller than its experimental value 92.4MeV, 
even if we use two-loop and assume very 
strong coupling in the IR region ao — 47r. This 
fact again suggests the importance of the scale 
ambiguity in the improved ladder SD equation. 

It should be emphasized that iVps = 0.8 is as- 
sumed irrespective to the value of ao in Figure |[ 
Since the factor A^ps is a decreasing function of 
ao, the ladder-exact value of f n is considered 
more stable than the result of the Pagels-Stokar 
formula. In order to make the analysis more re- 
liable, we therefore need to evaluate the ladder- 
exact value of /,r using the non-local gauge fixing 
method. Such a calculation is technically difficult 
to be performed, however. We thus leave the sub- 
ject as a problem to be examined in future. 
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FIG. 8: Results of (a) the pion decay constant and (b) the "constituent quark mass" m CO nst with use of 
the one- and two-loop MS couplings. The non-local gauge fixing method is adopted in the ladder SD equation. 
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FIG. 9: Results of (a) the pion decay constant and (b) the "constituent quark mass" m const with use 
of the effective coupling (two-loop RGE + finite part) with £ bg = 0. The non-local gauge fixing method is 
adopted in the ladder SD equation. The boundary condition of RGE is assumed to be af s (mb) = 0.2197, 
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FIG. 10: Results of (a) the pion decay constant f n and (b) the "constituent quark mass" m CO nst with use 
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We next turn to the case with the effective cou- 
pling, in which the scale ambiguity is expected 
to be resolved at the leading-order. Since we use 
the non-local gauge parameter £(z) Eq.(|o|) in 
our analysis of the SD equation, there does not 
exist a priori choice of the gauge parameter £bg 
in the effective coupling. We thus adopt both 
£b g = (Figure ||) and £ bg = 1 (Figure |l^) and 
compare the results. Again, the ao-dependence 
of the results is of non-negligible order. The cal- 
culated f v is relatively close to the observed value 
fir = 92.4MeV for a = 2-An with £ bg = 0, while 
£b g = 1 leads to a little bit larger predictions 
for the same range of ao- Considering that the 
non-local gauge £ (z) Eq.(|40|) approaches £ = in 
the asymptotic region, £b g = might be slightly 
better choice than £bg = 1- 



VII. SUMMARY AND DISCUSSIONS 

In this paper, we have calculated the pion de- 
cay constant f v from the high energy QCD cou- 
pling strength a ! f s (Mz) by using the improved 
ladder Schwinger-Dyson (SD) equation. The SD 
equation was analyzed both with and without 
the Higashijima-Miransky approximation for its 
angular integral. The non-local gauge param- 
eter method was adopted in the analysis with- 
out the Higashijima-Miransky approximation in 
order to keep the wave-function factor trivial 
A = 1. The effective coupling was calculated 
in the background gauge fixing method with ar- 
bitrary covariant gauge parameter £b g - Analyz- 
ing the Landau gauge improved ladder SD equa- 
tion combined with the £bg = next-to-leading- 
order effective coupling, we obtained / w = 85- 
106MeV depending on the value of af- s (M z ) = 
0.1172±0.0020 within the Higashijima-Miransky 
approximation. Our result impressively agrees 
with its experimental value / w = 92.4MeV and 
suggests quantitative validity of the improved 
ladder SD equation. It is interesting to compare 
our result on /„• with the previous analyses || |j 
where the value of the calculated is less than a 
half of its experimental value = 92.4MeV. In 
the previous analyses, the renormalization scale 
[i was naiv ely identified with the gluon momen- 
tum ^/]fc 2 | in the ladder SD equation. Such an 
identification has the problem of the scale am- 
biguity, however, as we pointed out in this pa- 
per. The improvement achieved in our analysis 
comes mainly from our use of the effective cou- 



pling, instead of the naive use of the one-loop 
MS coupling. The leading-order In fi dependence 
of a^ s (^) cancels with the ^-dependence of the 
finite correction in our effective coupling. 

We found that the effective coupling with 
£bg = 1 (the effective coupling derived in the 
pinch technique (PT)) leads to a larger value of 
Z^. This deviation of f n in £b g = 1 may be un- 
derstood as an indication of the double counting 
of the ladder diagrams in the analysis of the lad- 
der SD equation in £ = combined with the PT 
effective coupling (£b g = !)■ In the analysis of 
the non-local gauge parameter method, we were 
not able to obtain stable results. The results, 
however, were shown to be consistent with the 
results of the Higashijima-Miransky approxima- 
tion in order of magnitude. 

We next comment on a different approach to 
the scale ambiguity, in which a^- s (/i) is used with 
H tuned so as to minimize the finite correction. 
In the case of QED (or QCD in Np — > oo limit), 
such a scale is given by \i = -\/l^ 2 l exp(— 5/6) 
with k being the momentum of photon (gluon) 
propagator, while it is \i — i/[Pjexp(— 205/264) 
in the case of the gluon propagator (£ bg = 0) of 
Np = QCD. Since these scales are close to 
each other, it m ay b e possible to use universally 
the scale /i = ^/|fc 2 | exp(— 5/6) in the analysis of 
the improved ladder SD equation. It is straight- 
forward to evaluate the in such a simple pre- 
scription. Thanks to the scale invariant property 
of the SD equation, we just need to multiply the 
factor exp(5/6) to the result of /„• in the naive 
ladder SD equation with the MS coupling. Fig- 
ure |?] thus reads / w = 64-82MeV for the one-loop 
MS and /„ = 122-155MeV for the two-loop MS 
with this prescription. The difference between 
this method and the effective coupling method 
comes from the rest of the finite corrections and 
the treatment of the two-loop RGE. 

Many issues remain unsolved and need further 
investigation. In order to calculate the value of 
f n more precisely with the non-local gauge pa- 
rameter method, we need to evaluate the ladder- 
exact /„• by solving the Bethe-Salpeter equation 
of the pion. The problem of the £ b g dependence 
of the effective coupling should also be studied. 
It is rather non-trivial task to find adequate re- 
lation between £ (the gauge parameter in the 
SD equation) and £ bg (the gauge parameter in 
the effective coupling), especially in the non-local 
gauge parameter method. It is also interesting 
to investigate the critical behavior of the dense 
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and/or hot QCD using the effective coupling de- 
scribed in this paper. Since the method of the 
ladder SD equation is now smoothly connected 
with the high energy QCD, the high-density and 
high-temperature behaviors of the chiral phase 
transition can now be studied in a more trustful 
manner than before. 

Finally, the success of the ladder QCD implies 
that a bulk of driving force of the dynamical chi- 
ral symmetry breaking comes from the ladder- 
type diagrams. The result presented in this pa- 



per therefore provides deeper understanding of 
the low energy QCD dynamics. 
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